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ABSTRACT 

We determine the shape, multiplicity, size, and radial structure of superclus- 
ters in the ACDM concordance cosmology from z = to z = 2. Superclusters are 
defined as clusters of clusters in our large-scale cosmological simulation. We find 
that superclusters are triaxial in shape; many have flattened since early times 
to become nearly two-dimensional structures at present, with a small fraction of 
filamentary systems. The size and multiplicity functions are presented at differ- 
ent redshifts. Supercluster sizes extend to scales of ~ 100 — 200/i _1 Mpc. The 
supercluster multiplicity (richness) increases linearly with supercluster size. The 
density profile in superclusters is approximately isothermal (~ R~ 2 ) and steepens 
on larger scales. These results can be used as a new test of the current cosmology 
when compared with upcoming observations of large-scale surveys. 

Subject headings: large-scale structure of universe — cosmology: theory — galax- 
ies: clusters: general 
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INTRODUCTION 



Superclusters are the largest structures in the Universe. Their sizes span the range 
from a few times the radius of a typical galaxy cluster (~ 3-4 h~ x Mpc) to the lengths o f 
the "Great Walls" of galaxies obser ved by the CfA R edshift Survey (jGeller fc Huchralll989l ) 
and the Sloan Digital Sky Survey (IGott et al.ll2005l ). These latter sizes of ~ 10 2 /i _1 Mpc 
represent a significant fraction (~ 10%) of the horizon scale. 

Because they are so large relative to a typical survey volume, superclusters have proven 
difficult to study, especially in a statistical manner. Observationally, atte mpts to charac- 
terize the large-scale structure using distributions of galaxies date back to IHubbld (119361 ). 
but these metho ds did not be come powerful until the advent of galaxy cluster and galaxy 
redshift surveys. lAbelll (119581 ) used the Palomar Sky Survey to construct a catalog of rich 
clusters, with which he was able to find evidence for what he called "second-order clus- 



ters," i.e. clusters of c 



Gregory &: Thompson! ( 



uster s, which are indeed superclusters. Subsequent work included 
1978T)'s study o f the supercluster containing the rich Coma cluster, 



and lGregory. Thompson, fc Tifftl (1l98ll ) 's work on the Perseus supercluster. Other observa 



tional approaches to sup ercluster id e ntification have been based, like that o f 



on cluster catalogs (e.g. bortl Il983l: IBahcall fc Soneiral 11984 iBahcalll Il988l : 
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bj; lErdogdu et al. 



2004 
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19941: iKolokotronis. Basilakos. fc Plio nis 200 2]) or on analysis o f smoothed density fields 



Einasto et al.l 120061 ). These have been impor- 



tant in characterizing the low-redshift Universe, but they have not yet allowed the study of 
supercluster evolution from early cosmic times (i.e., from redshifts of order unity). How do 
superclusters form, and how do their shapes, multiplicities, sizes, and structures evolve with 
time? The answers to these fundamental questions are not yet known. 

Recently, there has been exciting progress o n both the observational and computatio nal 
fronts. The Sloan Digital Sky Survey (SDSS, lYork et al.l I2OO0I : Istoughton et al.l I2OO2I ) is 
collecting deep, broadband imaging and spectroscopy of a large, contiguous portion of the sky. 
SDSS has created the largest galaxy redshift survey cata log to date, and much work is also 



unde rway to improve photometric redshift techniques (e.g. lCsabai et all2000i : lBudavari et al. 



20001 ). which would permit measurements of 3-D galaxy positions for the much larger SDSS 
imaging catalog. On t he theoretical side , a concordance cosmological model has gaine d 
widespread acceptance (IBahcall et al.lll999l ; iBennett et al.ll2003l ; ISpergel et al.ll2003l . 120061 ). 
Applying this model to cosmological dynamics codes such as that described below, we can 
now reliably predict the properties and evolution of superclusters. The time is therefore 
ideal for a computational analysis of large-scale structure using new, powerful simulations, 
the results of which can be used to determine the properties and evolution of superclusters. 
In addition, these results can be compared with observations, thus providing a new test of 
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the current cosmological model — complementary to that provided by CMB, large-scale power 
spectra, and other studies — as well as shed light on structure formation and evolution. 

We begin by describing, in §2, the process by which clusters and in turn superclusters 
were identified in the simulation output. We provide quantitative descriptions of the mul- 
tiplicity and size distributions of the identified superclusters in §3.11 and §3.21 We examine 
the radial structure of our superclusters in §3.31 and the dimensionality of superclusters in 
§3.41 The evolution of the above properties with redshift is presented in £JD Results for 
superclusters selected from only higher-mass clusters are presented in SJ5J Our conclusions 
are summarized in Sj6l 



2. SELECTION OF SUPERCLUSTERS 



2.1. The Simulation 



We use the efficient, parallel dark matter simulation code of iBode fc Ostrikerl (120031 ). 
A large-scale, high-resolution simulation was evolved with a particle-mesh method to com- 
pute long-range gravitational forces, and a tree code to treat high-density regions. We 
use the cosmological parameters of the concor dance model: Q m = 0.27, fl^ = 0.73, H n 
70 km s" 1 Mpc" 1 , n s = 0.96, and as = 0.84 (cf. ISpergel et al.l (120031 ); see also ISpergel et al. 



(120061)). A total of 1260 3 particles, each assigned the mass 1.264 x 10 n /i _1 M o , is evolved 
in a periodic box lSOO/i^ 1 Mpc on a side, and positions are saved at the appropriate times 
in order to construct a light cone with vertex at one corner of the box. Clusters are then 
identified in two different, partially overlapping volumes: a "low- redshift" sphere of radius 
1500/i -1 Mpc, corresponding to a maximum redshift of z ~ 0.57, and a "high-redshift cone," 
which includes the low- z sphere's positive octant (x > 0, y > 0, z > 0) and extends this 
octant out to z = 3, a distance of ~ 4600/i _1 Mpc. Using the cluster selection scheme de- 
scribed below, a total of 1,442,616 clusters with M vir > 1.75 x lO 13 /i _1 M are found in the 
low- z sphere , and 720,550 in the high-^ cone. Further details on this light cone simulation 
are given by lHopkins. Bahcall. fc Bodd (120051 ). hereafter HBB. 



2.2. Cluster Selection 



Cluste rs are identified i n the simulation using the friends-of-friends (FOF) algorithm 
(see HBB; iDavis et al.lll985l ). The FOF results depend on the chosen linking length L. A 
common parameterization uses the linkage parameter b = L/r av , where r av is the mean 
interparticle separation. This is useful in applications where the mean density of objects, 
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and in turn their mean separation, changes over time. Holding b constant imposes the same 
minimum overdensity for clusters at all redshifts. We use the value b = 0.2 (see also HBB), 
corresponding to an overdensity of ~ 180 times the mean density. The center of each FOF 
cluster is taken to be the position of the most bound particle. This center is used to compute 
the virial radius R V i r , defined as the radius enclosing the virial overdensity expected from 
spherical top-hat collapse (which evolves with redshift in the ACDM cosmology). In this 
analysis, we consider only clusters with virial mass threshold (i.e., the total mass within the 
virial radius) M vir > 1.75 x 10 13 /i -1 M Q . This mass threshold is typical for a poor cluster of 
galaxies. The mean mass of all clusters in the low- z sphere is ~ 4.40 x 10 13 /i _1 M Q . In Sj5] 
we discuss superclusters selected from clusters with a higher mass threshold. 

The ellipticities of the clusters identified at low and high z are studied in detail by 
HBB. Relevant to our present work on superclusters, they find alignments between nearby 
cluster pairs out to separations of ^lOO/i" 1 Mpc; the strength of alignment increases with 
redshift and with decreasing cluster separation. HBB explain these alignments as being due 
to clusters forming and gathering material along "filamentary superstructures" ; they provide 
evidence for overdensities along the separation vectors between pairs of aligned clusters. 



2.3. Supercluster Selection 



As demonstrated by lAbelll (119581 ). iBahcalll (119881 ) and references therein, superclusters 



can be identified by their overdensity of clusters. We use the FOF method to identify super- 
clusters from the cluster sample defined above (§2.1-2.2). The mean density of clusters in the 
low- z sphere (we shall not discuss the high-z cone until §1]) is n av = 1.02 x lO -4 ^ -1 Mpc) -3 , 
yielding a mean intercluster separation of r av = riav^ = 21Ah~ 1 Mpc. We therefore select 
linking lengths of L = 3,5,7, and lOh" 1 Mpc, corresponding to b ~ 0.14,0.23,0.33, and 
0.47, respectively. All scales in this paper are in comoving coordinates. The different link- 
ing lengths sample different portions of the supercluster geometry. Small linking lengths 
identify only the dense cores of superclusters, whereas larger linking lengths allow the FOF 
algorithm to percolate to lower-density filaments of clusters. Large L-values therefore reach 
to the outer parts of superclusters, as well as identify looser superclusters, and often connect 
structures that appear discrete with smaller L. 

We use Monte Carlo simulations to compare the significance of the selected superclusters 
with that of random associations. Clusters are placed at randomly chosen positions within 
a sphere such that the total number density is the same as that for the simulated clusters, 
and the FOF method is applied, as before, with L = 3, 5, 7, and lOh^ 1 Mpc. We find 
significantly fewer clusters per supercluster, and significantly fewer superclusters down to 
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multiplicity as small as 3, for all four linking lengths tested. This result implies that most of 
the superclusters are real, even for L = lOh^ 1 Mpc. The only exception is the case of binary 
clusters (superclusters with only 2 members) selected with L = \0ti~ 1 Mpc, of which we find 
only ~ 25% more using the simulated clusters than we do with the Monte Carlo sample. 
The random rate of superclusters is typically < 10% for iV > 3 superclusters, decreasing 
to < 1% for N > 6 superclusters. The Monte Carlo results are compared with the ACDM 
multiplicity functions below. 

To help visualize the typical geometries of superclusters identified by the FOF method, 
we present supercluster maps of a typical region (a cube of side length 200/i~ 1 Mpc, concentric 
with the low- z sphere) for L = lhr x Mpc (Figs. [012]). Each point is a cluster, and a given 
symbol represents membership in the same supercluster. The maps present a projection of 
the 3-D distribution onto the xy-plane in the first plot, and onto the xz-plane in the second 
plot. Viewing the plots in turn corresponds to a 90° rotation of the viewing angle about the 
x-axis. For clarity, we show only superclusters with at least 5 member clusters. Members of 
a supercluster that extend outside the limits of our cubical region are not shown, so some 
superclusters appear to have fewer members than the imposed minimum. Comparing these 
maps to analogous figures generated with L = 5 and lOh^ 1 Mpc reveals, as expected, that 
increasing the linking length often causes neighboring superclusters to be connected into 
one. This is especially true when L increases from 7 to lOh^ 1 Mpc, resulting in many large, 
complicated structures with central cores joining multiple lower-density filaments. 



3. SUPERCLUSTER PROPERTIES 

3.1. Multiplicity Function 

The multiplicity of a supercluster is the number of clusters it contains, N. In Fig.[3]we 
present, for various linking lengths, the integrated supercluster multiplicity function — the 
number density (per cubic comoving h~ x Mpc) of superclusters with N or more members, 
as a function of N. Also shown for comparison are the corresponding results from random 
catalogs, as described in the previous section. 

Longer linking lengths yield more superclusters of all multiplicities, as expected. Shorter 
linking lengths yield curves that rise more steeply toward low multiplicities, implying that 
they preferentially find small superclusters with low multiplicity. Conversely, longer linking 
lengths find a higher percentage of larger, richer superclusters. The shorter linking lengths 
naturally trace high-density regions (e.g., dense cores of superclusters), while larger linking 
lengths characterize lower-density regions, such as the outer parts of superclusters or dif- 
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fuse superclusters. Small L can also divide structures that are identified with larger linking 
lengths. The fact that we see superclusters with ~ 40 or more clusters may seem surpris- 
ing, but the total number density of such structures is only 10 -8 per {h~ l Mpc) 3 , and it is 
based on a very low mass threshold of cluster members, M vir > 1.75 x 1O 13 /i~ 1 M . Such 
clusters are numerous, and any comparison with observations should of course account for 
this threshold. For superclusters selected using higher-mass clusters, see §4.11 We repeated 
the above analysis after selecting superclusters in redshift-space coordinates; the resulting 
multiplicity functions are nearly identical to the real-space functions at low multiplicities, 
and are shifted to slightly higher N (by up to 20%) at the highest multiplicities. 



3.2. Size Function 



The supercluster size function — i.e., the number of superclusters above a given size - 
is presented in Fig. H] for different linking lengths. The size of a supercluster is defined as the 
maximum distance of any member cluster from the supercluster center of mass. As with the 
multiplicity function, the size function rises more steeply for shorter linking lengths. This 
implies that the longer linking lengths find a greater ratio of large superclusters to small 
superclusters, analogous to the result for the multiplicity function. Similar results are found 
using redshift-space coordinates; for the largest linking length used, the sizes are slightly 
(< 5%) larger. 

We find that the largest supercluster radii are ~ 80/i _1 Mpc, corresponding to diameters 
of ~ 160/t" 1 Mpc. This is consistent with the l argest structures seen in observations (e.g. 
Bahcall fc Soneiral Il984j ; iGeller fc Huchral Il989l ; iGott et al.l 120051 ) , and it is also the largest 
distance at which HBB detect cluster alignment in the simulation. 



3.3. Radial Structure 

The average multiplicity of superclusters is plotted as a function of supercluster size for 
different linking lengths in Fig.|5J For N = 2 or 3 — objects more aptly called binary or triple 
clusters than superclusters — the maximum radius is set mainly by the separation(s) of the 
2 or 3 clusters; thus the average multiplicity changes slowly over this range of separations. 
Beyond this regime, the multiplicity and size of superclusters appear to be related nearly 
linearly. This result suggests that the spherically averaged density profile in superclusters is 
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roughly isothermal] i.e., n(r) oc r 2 , where the cluster number density n(r) satisfies: 



This profile yields multiplicity N oc R max for n(r) oc r 2 . 

We can investigate the spherically averaged density profile more directly by examin- 
ing the distribution of clusters within individual superclusters and determining the average 
number density profile for the entire sample. To do so, we superpose the centers-of-mass 
(calculated using only the mass in clusters) of all superclusters, and scale all superclusters 
to the same size by dividing all cluster coordinates (relative to the center-of-mass) for a 
given supercluster by its R max - We add up the total number of clusters inside radial bins, 
divided by the total number of clusters in superclusters for the given linking length. We have 
included only superclusters with 5 or more members; because the multiplicity function drops 
off rapidly, most of the superclusters included have only 5 or 6 members. This fact affects the 
combined profile in two ways. First, only in few cases will one of the five members happen 
to lie very near the supercluster's center-of-mass; thus it is only at R/R max > 0.2 that we 
begin to find a significant fraction of the clusters. Second, one out of roughly five clusters in 
each supercluster lies at R/R max = 1- Thus, for the smaller linking lengths, the integrated 
number fraction approaches a limit of ~ 0.8; for the larger linking lengths, which find richer 
superclusters, the limiting fraction is slightly higher. The fraction of clusters within R/R max 
increases linearly with radius for the range 0.3 < R/R max < 0.7. This linear relationship 
indeed implies n(R) ~ R~ 2 for the supercluster density profile in this range. The profile 
steepens on larger scales. 

The average mass of clusters in superclusters, M = 4.8 x lO 13 ^" 1 M Q for the L = 
7hr x Mpc superclusters, is somewhat larger than the average mass of clusters in the total 
sample, M = 4.4 x 10 13 /i _1 M , indicating that superclusters contain slightly more massive 
clusters, on average, than the mean. The L = 7/i~ 1 Mpc superclusters with N > 5 members 
have an even higher average cluster mass of M — 5.5 x 10 13 /i _1 M , implying that richer 
superclusters contain more massive clusters than poorer superclusters. 



We measure the shape, and in particular the dimensionality, of a given supercluster by 
fitting a 3-D ellipsoid to the distribution of its member clusters. We use the code developed 
by HBB to fit ellipsoids to the dark matter particles comprising their clusters. The code 
constructs, for a given supercluster, the 3x3 matrix of second moments of cluster positions 




(1) 



3.4. Supercluster Shapes 
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relative to the supercluster center-of-mass, i.e. 

hj = XiXjin, (2) 

where Xj and Xj are two of the position coordinates for a given cluster, and m is the mass 
of the cluster; the sum is over all member clusters. The eigenvalues of this matrix are 
simply the three axis lengths of the best-fit ellipsoid for the supercluster, times a known 
constant. Following the notation used by HBB, we shall denote these lengths 01,02,03 for 
the primary, secondary, tertiary axes, respectively; i.e., a\ > a 2 > 03. We restrict our 
analysis to superclusters with N > 5; for poorer superclusters, the "best-fit ellipsoid" would 
not be very meaningful. 

We have already examined size distributions for superclusters ( §3.2[) . so we now focus 
on the axis ratios, which allow us to probe the shape and dimensionality of superclusters. 
Fig. O displays the distribution of primary axis ratios a 2 /ai for superclusters selected with 
L = 5, 7, 10o _1 Mpc. Fig. [7] shows the distribution of secondary axis ratios o 3 /a 2 , and Fig.[H] 
shows 03/01. The mean and peak axis ratios from each curve in these figures are listed in 
Table 1. 

The secondary axis ratios are generally smaller than the primary axis ratios. We find 
that a primary axis ratio of a 2 /oi « 0.6 is typical for superclusters selected using a broad 
range of linking lengths; this corresponds to a primary ellipticity of ei = 1 — a 2 /a\ ~ 0.4. 
The secondary axis ratio is more strongly dependent on the linking length used to find 
the superclusters. For L = 5/i _1 Mpc, we find a peak 03/02 of ~ 0.4, corresponding to 
a secondary ellipticity of e 2 = 1 — 03/02 ~ 0.6; for L = lO/T^Mpc, we find a peak of 
03/02 ~ 0.2, or t 2 ~ 0.8. The peak 03/01 ratios range from « 0.13 for L = 10/i _1 Mpc to 
« 0.23 for L = 5h~ l Mpc. 

These measurements allow us to construct a picture of the typical shape and dimen- 
sionality of superclusters. The interpretation of different ellipticities can be described as 
follows: 

1 . A combination of high primary and low secondary axis ratios is typical of two-dimensional 
pancake-like structures (two dimensions are large, and one is small). 

2. A combination of low primary and high secondary axis ratios is typical of one- dimensional 
filamentary structures (one dimension is large, and two are small). 

3. A combination of primary and secondary axis ratios that are both significantly less than 
unity is typical of structures most appropriately described as triaxial (each dimension 
is a different size). 
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4. If both axis ratios are close to unity, then the distribution is nearly spherical. 

The results we find suggest that superclusters are typically triaxial structures nearing 2- 
D, pancake-like structures at z ~ 0. They are neither spherical nor purely filamentary 
in nature. As the linking length is increased, the resulting superclusters become more two- 
dimensional; the cores of superclusters are thus more triaxial, spreading out to more pancake- 
like structures on larger scales. 

We plot the bivariate distribution of a^/ai and ct 3 /a 2 for L = llr 1 Mpc in Fig. [9] (the 
distributions for L = 10 and 5/i" 1 Mpc are similar). Applying the interpretation rules listed 
above, 1-D structures reside in the upper left corner of Fig.El 2-D in the lower right corner, 
and 3-D (spherical) in the upper right. Points in the middle or toward the lower left corner 
represent triaxial structures. 

Fig. IH] shows that a large majority of the superclusters are more nearly two-dimensional 
than one-dimensional. No nearly spherical superclusters are found with this linking length, 
but there is a smaller population of filamentary structures. Similar distributions are found 
for L = 10 and 5o _1 Mpc, though showing more flattened structures (2D) as L increases, 
and more triaxial structures as L decreases. This is the same trend depicted in Fig. [7] as a 
shift toward lower secondary axis ratio with increasing linking length. The shorter linking 
lengths also yield significantly higher fraction of filamentary structures — i.e., superclusters 
with low primary axis ratio and a secondary axis ratio approaching unity. 

This trend suggests that most superclusters at z ~ are triaxial structures with one 
dimension considerably smaller than the other two. This geometry is apparent in Fig. El and 
becomes more significant when a larger FOF linking length is used. Since the superclusters 
are not completely flat, smaller linking lengths may probe only the dense core regions, which 
are smaller in size than the supercluster thickness. Thus, while the thickness sets an upper 
bound to the tertiary axis length identified with any linking length, the higher L-values find 
longer primary and secondary axes by probing farther out into the pancake. This causes 
larger linking lengths to find a smaller mean secondary axis ratio, and reveal the underlying 
2-D nature of the large-scale superclusters. Similarly, for a given L, the largest, highest- 
multiplicity superclusters (N > 10) are on average flatter (03/02 ~ 0.24 for L = \0h~~ 1 Mpc) 
than those with N = 5-10 (03/02 ~ 0.35). This flatness is still apparent in projection: 
for superclusters selected (from M vir > 1.75 x lO 13 /^ 1 M clusters) in redshift-space, the 
projected axis ratios have a broad distribution with a mean 02/01 = 0.4 ± 0.005 and a peak 
02/01 ~ 0.35, for all linking lengths L =5, 7, and lOo" 1 Mpc. 

The triaxiality of superclusters can also be quantified using the measurement commonly 
employed for geometrical descriptions of elliptical galaxies, first introduced explicitly by 
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Statlerl (J199J): 

= 1 ~ (Q2/Q1) 2 

" 1 - (a 3 /ai) 2 1 J 

With a\ > ci2 > ct3, T will take on values between and 1 (unless ai = a 2 = «3j in which 
case T is undefined). Two-dimensional pancake structures approach T — ► 0, while filaments 
approach T — > 1 (note, however, that the increase is not linear with axis ratio). Intermediate 
values of T represent triaxial systems, with triaxiality increasing with T. By this measure 
our low-z superclusters are triaxial, ranging from a mean of T fa 0.65 for L = lOh" 1 Mpc to 
T « 0.69 for L = Mpc. 

For all linking lengths we see a tail in the distribution that describes structures which 
are very nearly one- dimensional. Thus, some mass concentrations are filamentary at z ~ 
0, especially in the relatively dense regions probed by L = bh" 1 Mpc. Indeed, we would 
expect structures that have collapsed along two dimensions to be of higher average density 
than those that have collapsed along only one dimension. This result agrees with HBB's 
identification of filamentary density enhancements in the large-scale structure; here we add 
the prediction of more diffuse flattened structures surrounding and possibly connecting the 
filaments. 



4. SUPERCLUSTER EVOLUTION 

We examine the time evolution of the supercluster properties studied above. We study 
superclusters in two redshift slices at 0.8 < z < 1.2 and 1.7 < z < 2.4, denoted z ~ 1 and 
z ~ 2, respectively. 

Superclusters are identified by the FOF method as before. However, since the cluster 
abundance decreases at higher redshift, and their mean separation thus increases, we compare 
superclusters found at different redshifts not with the same linking length L, but instead 
with the same linkage parameter b = L/r av , where r av is the mean co moving intercluster 
separation. As discussed in §2.2, this approach sets a minimum overdensity for detection of 
superclusters, which is constant with redshift. 

This use of a constant linkage parameter accounts for the reduced density of clusters at 
high redshift. However, because clusters are clustered, the density of superclusters depends 
on the relative rates of cluster formation in regions of high density vs. low density. Hierar- 
chical supercluster formation, in which clusters gravitate toward each other over time, also 
contributes to a reduced supercluster density at early times. 

We indeed find a lower density of superclusters in the high-z cone. The total number 
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of superclusters decreases by a factor of a few from z ~ to z ~ 1, and by over an order 
of magnitude from z ~ 1 to z ~ 2. In order to construct a sufficiently large sample, we 
consider broad ranges in redshift, especially for the z ~ 2 slice. In addition, we omit the 
6-value corresponding to the L = 3k' 1 Mpc (z fa 0) case from our subsequent analysis, as it 
results in too few superclusters at high redshift. 

As stated in §2.3, our z ~ linking lengths of L = 5, 7, 10/i _1 Mpc correspond to 
b = 0.23, 0.33, 0.47, respectively. We use these latter linkage parameters b to refer to the 
linking lengths used in our subsequent analysis. A given value of the 6-parameter can be 
converted to an actual comoving linking length at a given redshift simply by multiplying by 
the mean intercluster separation at that redshift. The mean cluster densities in our z ~ 1 
and z ~ 2 shells are n av = 3.53 x 10 -5 and 2.57 x 10 _6 (/i _1 Mpc)~ 3 , corresponding to mean 
cluster separations of r av = 30.5 and 73.0k' 1 Mpc, respectively. Therefore, b = 0.23, 0.33, 
and 0.47 correspond to L = 7, 10, and 14.3/i" 1 Mpc at z m 1, and 16.8, 24, and 34.3/i- 1 Mpc 
at z ~ 2. 



4.1. Multiplicity Function 

The evolution of the multiplicity function is presented in Fig. [TO! We present the in- 
tegrated multiplicity function for all three redshifts, using two different values of b. The 
number density of superclusters drops noticeably with increasing redshift. Monte Carlo sim- 
ulations analogous to those described in §2.3 were also run at z « 1 and 2, with results 
similar to those obtained at low z: the multiplicity functions are physically significant for 
all multiplicities and linking lengths, except for binary clusters in the b = 0.47 case, where 
only ~ 30% more binary clusters are found over the Monte Carlo results. 

The redshift evolution of the multiplicity function depends on the linking length of the 
superclusters. The vertical distances between the b = 0.47 curves in Fig. [10] increase with 
increasing multiplicity, implying that the abundance of the richest (highest-iV) superclusters 
found with this linking length drops off especially swiftly with increasing z. However, for 
b = 0.23 the vertical distances decrease with increasing N, to the extent that the abundance 
of the richest superclusters found with this linking length is in fact roughly the same at z « 1 
as it is today. We interpret these trends as further evidence that small linking lengths (e.g. 
b = 0.23) select the densest regions of superclusters. Whereas the larger linking lengths do 
not find many rich superclusters at high redshift because of the lack of clusters, the smaller 
linking lengths find the same dense supercluster cores even at much earlier times. These 
dense cores are likely the sites of the largest overdensities present at early times, and the 
first places where observable structures form through gravitational collapse. Thus smaller 
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linking lengths relative to the mean intercluster separation are needed to locate the highest 
overdensities in the observable Universe. The lack of rich superclusters at high z is also 
evidence that these structures form largely through gravitational accretion of clusters that 
are initially too distant to be identified as supercluster members. 

Supercluster size remains highly correlated with multiplicity at high z, but with a lower 
iV — R max normalization than at low z. The size function does not evolve significantly with 
redshift for the larger superclusters (see Fig.[TT]), while, as expected, the number of cluster 
members decreases. 



4.2. Dimensionality 

Using the method described in §3.4[ we fit ellipsoids to the superclusters identified 
at z ~ 1 and z ~ 2. Because we only perform these fits on superclusters with iV > 5, 
which are rare at high redshift, the statistical uncertainties are larger. For example, only 68 
superclusters are identified at z ~ 2 for b = 0.23. 

We use the axis ratios of the best-fit ellipsoids to characterize the evolution of super- 
cluster dimensionality. Tables 2 and 3 present the mean primary and secondary axis ratios, 
respectively, for b = 0.23,0.33,0.47 at z ~ 0,1,2. We also construct the axis ratio distri- 
butions (similar to Figs.EHT]) as a function of z for b = 0.33 and 0.47 (there are too few 
b = 0.23 superclusters at high redshift to create meaningful distributions). The results for 
b = 0.47 are shown in Figs. [T2irT3l The number of high-redshift superclusters is too low to 
construct bivariate distributions as in Fig. [HI but we can consider how the changes in mean 
ratios correspond to translations of these distributions. 

The strongest evolution observed is the increase in the secondary axis ratio with in- 
creasing redshift, for b = 0.47 (Fig. [T3|) . Most of the change occurs between z ~ and z ~ 1, 
with only a slight additional increase between z ~ 1 and z ~ 2 (see Table 3). Typical sec- 
ondary axis ratios increase from a 3 /a 2 ~ 0.25 at low 2: to a 3 /a 2 ~ 0.45 at high z. Thus the 
flattened structures observed with b = 0.47 at z ~ (see §3.41) appear less flattened at high 
redshifts, indicating their gravitational collapse over time. Very little evolution is evident 
in the secondary axis ratio for superclusters identified with b = 0.33, but the population of 
high-density filaments with a 3 /a 2 ~ 1 evolves between z « 1 and z ~ 0, as there are almost 
none of these filamentary structures at high redshift. The z = distributions for b = 0.33 
are shown in Figs.lMZl & t higher redshifts the distributions are similar to the corresponding 
curves in Figs. [T2lfT3l 

The primary axis ratios are slightly smaller at higher redshifts for both b = 0.33 and 
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0.47. Typical values decrease from 02/01 ~ 0.6 at z ~ to 02/01 ~ 0.5 at z ~ 1-2. In a 
plot similar to Fig.[9], this would correspond to a leftward shift in the distribution, i.e. a 
shift toward greater triaxiality at high redshift, and a shift to a 2-D geometry as the systems 
approach z ~ 0. Quantitatively, our high- z superclusters have mean triaxiality T « 0.79 for 
both linking lengths, significantly higher than the mean T for any linking length at low z. 

We conclude that superclusters had highly triaxial structure at early times, and that 
between a cosmological look-back time of ~ 8 Gyr (corresponding to z ~ 1 for the cosmology 
described in §2.1) and the present day, many collapsed along their smallest dimensions to 
form the more nearly pancake-like structures observed in the simulation at z ~ 0. In the 
same span of time, some structures in regions of highest density also collapsed along a second 
dimension to form nearly one-dimensional filaments. 



5. SUPERCLUSTERS SELECTED USING HIGH-MASS CLUSTERS 

How do the supercluster properties change if superclusters are selected using higher- 
mass clusters as their seeds? While the abundance of high-mass clusters is significantly 
lower, and the statistical uncertainties larger, the massive clusters are more easily detected 
in observational surveys. We repeat the analysis described above for superclusters selected 
using the rarer, higher-mass cl usters with M,„> > l xlO 1 4 ^" 1 M . This threshold corresponds 
to typical "rich clusters" (e.g. lAbelll[l958l : lBahcaJ]|ll988h . There are 9.68 x 10 4 clusters above 



this mass threshold at z ~ 0, and 4.88 x 10 3 clusters at z w 1. The cluster abundance 
is 6.8 x lO-^/T^Mpc)- 3 at z « 0, decreasing to 0.8 x lO" 6 ^ 1 Mpc)~ 3 at z « 1. The 
intercluster mean separation is 52.7h~ l Mpc at z ~ 0, increasing to 107.9o _1 Mpc at z « 1. 

Superclusters are selected from these clusters at z ~ and z ~ 1 using the linkage 
parameters b = 0.33 and 0.47; these correspond to L — 17A and 24. 8/1 -1 Mpc respectively at 
z ~ 0, and L = 35.6 and 50.7o _1 Mpc at z ~ 1. The large linking lengths L reflect the large 
mean intercluster separation of these massive clusters. We also include, for comparison, the 
smaller linking length of L = 10o _1 Mpc at z ~ 0. There are fewer superclusters identified, 
and a smaller number of cluster members in each, as compared to those selected using the 
low-mass clusters. The overall structure of the superclusters is not significantly changed, as 
described below. 

The multiplicity function of the superclusters is presented in Fig. [14] for both z ~ 
and 1. The results show, as expected, a similar shape to the multiplicity function of su- 
perclusters selected using the more numerous low-mass clusters (Figs.O [101 [TT|) . but with 
a greatly reduced amplitude: the number of superclusters at a given multiplicity is lower, 
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and — equivalently — the supercluster multiplicity (richness) is smaller for a given supercluster 
abundance. This is partly due to the much smaller number of massive clusters. The evolu- 
tion of the multiplicity function is stronger than that of superclusters selected with low-mass 
clusters (cf. Fig. [TO]). Superclusters with linkage parameters of b ~ 0.2 — 0.33 contain up to 
~ 10 — 15 rich clusters at z ~ 0. The supercluster multiplicity decreases sharply at higher 
redshifts. 

The supercluster size function is presented in Fig.[Tj3 The maximum size of z ~ super- 
clusters ranges from ~ 40o _1 Mpc (i.e., twice the plotted R max radius) for L = 10o _1 Mpc, 
to ~ 90A _1 Mpc for L = 17 Ah- 1 Mpc {b = 0.33), to ~ 200k- 1 Mpc for L = 24.8o~ 1 Mpc 
(b = 0.47). For a given linkage parameter b, the maximum supercluster size is somewhat 
larger than that of superclusters selected with low-mass clusters, since the linking length 
values L are larger. For a given value of b, no significant evolution is observed in the size of 
the largest superclusters, reflecting the existence of these extended structures at early times. 
The same result is observed for superclusters selected with the lower mass threshold clusters. 
The number of cluster members within these large structures decreases with redshift, as seen 
above, due to the lowered cluster abundance at high z. The supercluster multiplicity versus 
size relation is similar — but reduced in amplitude, as expected — to the superclusters selected 
with the lower mass threshold (Fig.EJ). A nearly linear relation between multiplicity and size 
is observed for all but the binary clusters. 

The ellipticity distribution shows that superclusters exhibit a triaxial geometry on av- 
erage, but with significantly larger scatter than the more numerous superclusters selected 
with lower-mass clusters (Figs. fTSlfTB"]) . The mean primary and secondary axis ratios are each 
~ 0.4 to 0.5; for example, for b = 0.33 at z ~ 0, the mean values are 02/01 = 0.46 ±0.01 and 
a ?J&2 — 0.40 ± 0.01. These ratios do not change significantly with either linking length or 
with redshift to z » 1. The distribution of ellipticities around this mean is, however, quite 
broad. The supercluster shapes are not significantly different from those selected using the 
lower-mass clusters, except for the stronger flattening observed in the latter at z ~ 0, where 
the mean triaxial shape has evolved slightly towards a flattened two-dimensional shape (with 
a 2 /ai ~ 0.6 and a 3 /a 2 ~ 0.3 — 0.4). The superclusters selected with high-mass clusters are 
triaxial, with somewhat higher triaxiality at z ~ than those selected using lower-mass 
clusters — i.e., the high- mass superclusters are approaching a slightly more filamentary na- 
ture, with mean axis ratios of ~ 0.4 in both 02/01 and 03/02 (cf. Fig.[9]), and a mean 
triaxiality of T = 0.8 - 0.85. 
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6. CONCLUSIONS 

We have investigated the multiplicities, sizes, radial structures, and dimensionalities of 
superclusters from z ~ to z ~ 2 in ACDM cosmology. The results improve our under- 
standing of the properties (especially shapes) and evolution of superclusters, and can also 
be used as predictions for testing the cosmology with future observations. 

At low redshift, we predict a total number density of ~ 4 x 10~ 6 , 2 x 10~ 6 , and 
3 x l§~ 7 {h~ l Mpc) -3 superclusters with at least five member clusters with mass M v i r > 
1.75 x 10 13 /i _1 M for superclusters selected with linking length L = 10,7, and 5/i _1 Mpc, 
respectively. The maximum size of these superclusters ranges from ~ 150/i _1 Mpc for 
L = 10/i -1 Mpc to ~ 30/i -1 Mpc for L = 5/i _1 Mpc. The abundance of superclusters de- 
creases rapidly with increasing redshift. There are many more superclusters at low z than 
at high z largely because there are many more clusters that have formed by the present 
time and that have gravitated to form superclusters. The richest superclusters are almost 
entirely missing at high redshifts; these structures form late in cosmic time. However, the 
densest "core" regions of superclusters, traced well by FOF linking lengths ~ 25% of the 
mean intercluster separation, are present at z « 1 almost as abundantly as they are today. 
These are likely the sites of the highest overdensities in the early Universe. 

The spherically averaged density profiles of superclusters are well fit by an isothermal 
profile, n(r) ~ r~ 2 , over a broad range in radius. The profile is shallower at small radii 
and is steeper at large radii. A nearly linear relation exists between supercluster size and 
multiplicity, arising from this density profile. 

The clusters residing in superclusters are more massive on average than un-clustered 
clusters, and clusters residing in rich superclusters are more massive than those in poorer 
superclusters. 

We find that superclusters are triaxial in shape, especially at early times, where the mean 
triaxiality of our sample is T « 0.79. Over time, many collapse along a single dimension to 
approach two-dimensional shapes on the largest scales (i.e., those sampled by the FOF linking 
length L = \01i~ 1 Mpc at z ~ 0, using low-mass clusters with M„ r > 1.75 x 10 13 /i _1 M ); 
their cores remain triaxial. Quantitatively, the mean triaxiality for long linking lengths 
decreases to T w 0.65 at z ~ 0, whereas for smaller L (i.e., in the cores) it decreases to 
T « 0.69. The mean primary axis ratio of superclusters at z ~ is ~ 0.6, corresponding 
to a primary ellipticity t\ ~ 0.4. Typical secondary axis ratios at z 9a are 02/0-1 ~ 
0.45 and 0.3/0,1 rj 0.25 for L = 5/i _1 Mpc, and 02/01 ~ 0.25 and 03/01 ~ 0.12 for L = 
10/i _1 Mpc. Gravitational collapse increas es ellipticity, as the pote ntial gradient is larger 



along the minor rather than the major axis (ILin. Mestel. fc Shulll965l ); collapse thus leads to 



-16- 



lower dimensionality (cf. lYoshisato et al.ll2006l . for a more recent treatment). A nonnegligible 
population of dense, filamentary superstructures is also present at z ~ 0, though not at higher 
redshift. 

Superclusters selected using higher-mass clusters, M vir > 10 14 /i _1 M , show consistent 
results, with a reduced multiplicity function, as expected. These superclusters are triaxial 
in shape, with somewhat greater mean triaxiality at z ~ (T fh 0.8-0.85) than superclusters 
with lower-mass clusters. 

These properties, derived from a large-scale, high-resolution cosmological simulation, 
provide direct information on the size, content, shape, and evolution of superclusters, as 
well as testable predictions of the currently favored ACDM model with which future and 
currently accumulating observations can be compared. Once SDSS is complete, the best sta- 
tistical analyses to date of large-scale structure using observations will be possible, especially 
if accurate and well-behaved photo- z algorithms are develope d. There are other exciting ob- 
serving programs on the horizon. The Pan-STARRS Project (IKaiser et al.ll2002l ) is expected 
to see first light in 2006, and will carry out multi-band imaging surveys to identify galaxy 
clusters, as well as gravitational weak lensing observa tions to map the matter distribution in 
the Universe. The Large Synoptic Survey Telescope (jTyson 2002 ). which is expected to see 
first light in 2012, will provide improved weak lensing maps that will probe the dark matter 
distribution to the largest scales. The analysis presented here will aid in the comparison of 
these major observational data with the predictions of the concordance cosmological model. 
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Table 1: Supercluster Axis Ratios at z ps 0, Cluster Masses M vir > 1.75 x 10 13 o 1 M, 



L 


Mean 02/01 


Peak 02/01 


Mean 03/02 


Peak 03/02 


Mean 03/01 


Peak 03/01 


5 


0.60 ± .009 


0.64 


0.50 ± .007 


0.41 


0.29 ± .004 


0.23 


7 


0.57 ± .004 


0.62 


0.43 ± .003 


0.31 


0.23 ± .001 


0.19 


10 


0.54 ± .002 


0.64 


0.33 ± .001 


0.20 


0.16 ± .001 


0.13 



Note. — Peaks and mean values of the distributions shown in Figs.EHHl Statistical errors are given for the 
mean values. 



Table 2: Evolution of the Primary Axis Ratio a 2 /oi, M vir > 1.75 x lO 13 /^ 1 M, 



z 


b = 0.23 


b = 0.33 


b = 0.47 





0.60 ± .009 


0.57 ± .004 


0.54 ± .002 


1 


0.44 ± .03 


0.47 ± .01 


0.48 ± .01 


2 


0.52 ± .06 


0.50 ± .03 


0.50 ± .02 



Note. — Mean values are given as a function of the linking length parameter b. Statistical errors are 
given. Note that b = 0.23 corresponds to L z= o = bh^ 1 Mpc, b — 0.33 to Lq — Th^ 1 Mpc, and b — 0.47 to 
Lo = 10/i- x Mpc. 



Table 3: Evolution of the Secondary Axis Ratio 03/02, M V i r > 1.75 x 10 13 o 1 M 



z 


b = 0.23 


b = 0.33 


b = 0.47 





0.50 ± .007 


0.43 ± .003 


0.33 ± .001 


1 


0.36 ± .02 


0.40 ± .01 


0.44 ± .01 


2 


0.41 ± .05 


0.41 ± .03 


0.46 ± .01 



Note. — Mean values are given as a function of the linking length parameter b. Statistical errors are 
given. 
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X-Y Projection of z = Superclusters with L = 7 Mpc/h, N >= 5 




x 



Fig. 1. — Clusters identified as supercluster members at redshift zero, for L = 7h~ l Mpc and 
minimum multiplicity 5. Clusters represented by similar points are members of the same 
supercluster. See the text (§2.3) for further description. 



X-Z Projection of z = Superclusters with L = 7 Mpc/h, N >= 5 




-100 1 1 1 1 1 

-100 -50 50 100 



x 

Fig. 2. — Same parameters as Fig.fU but a different projection. 




Fig. 3. — Integrated supercluster multiplicity functions for linking lengths L = 3,5, 7, and 
10/i _1 Mpc, all at z ~ 0. Curves show superclusters identified in our random Monte Carlo 
simulations to test the physical reality of the superclusters. 
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Fig. 4. — Integrated size functions for linking lengths L = 3, 5, 7, and 10/i - 1 Mpc at z ~ 0. 
The size measurement along the horizontal axis is, for a given supercluster, the distance of 
the member cluster most distant from the supercluster center of mass. All are comoving 
scales. 



-23- 



100 - 



Supercluster Multiplicity vs. Size 
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Fig. 5. — Multiplicity vs. size of superclusters, z ~ 0. -R max is defined as in Fig.Hl Linking 
lengths L are quoted in Mpc. 
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Fig. 6. — Distributions of supercluster primary axis ratios for various linking lengths, at 
z « 0. Linking lengths L are quoted in Mpc. 
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Fig. 7. — Distributions of supercluster secondary axis ratios for various linking lengths, at 
z ~ 0. Linking lengths L are quoted in h~ l Mpc. 
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Fig. 8. — Distributions of supercluster third axis ratios 03/01 for various linking lengths, at 
z ~ 0. Linking lengths L are quoted in h~ l Mpc. 




Fig. 9. — Bivariate distribution of primary and secondary axis ratios for low- z superclusters 
found with linking length L = 7h~ l Mpc. See the text ( §3.4[) for interpretation. 



-28- 



0.0001 f 



o 



10- 



10- 



io- 



s io- B - 



A 



10" 



10 -io _ 



10" 



n 

o 



o 



Supercluster Multiplicity Function 



O 



o 



M cl ^1.75xlO 13 h- 1 M 

b = 0.47 

z = (L=10) • 
z = l (L=14) 
z = 2 (L=34) + 

b = 0.23 

z = (L=5) 
z=l (L=7) 
z = 2 (L=17) O 



O 








o 



10 



N 



100 



cl/SC 



Fig. 10. — Evolution of the integrated supercluster multiplicity function with redshift. Link- 
ing lengths L are quoted in h' 1 Mpc. 
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Fig. 11. — Evolution of the integrated supercluster size function (comoving scales) with 
redshift. Linking lengths L are quoted in h~ l Mpc. 
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Fig. 12. — Evolution of supercluster primary axis ratios with redshift. All curves are gen- 
erated from superclusters selected with linking length parameter b = 0.47. The points are 
connected for clarity. Linking lengths L are quoted in h^ 1 Mpc. 
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Fig. 13. — Evolution of supercluster secondary axis ratios with redshift. All curves are 
generated from superclusters selected with linking length parameter b = 0.47. Linking 
lengths L are quoted in h^ 1 Mpc. 
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Fig. 14. — Integrated supercluster multiplicity functions for superclusters selected using only 
high-mass clusters, M vir > 10 14 /i -1 M Q . Linking lengths L are quoted in h~ l Mpc. 
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Fig. 15. — Integrated supercluster size functions for superclusters selected using only high- 
mass clusters, M vir > lO 14 ^" 1 M . Linking lengths L are quoted in h~ x Mpc. 



